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Preface 


This is the third quarterly report on Contract NAS5-21798, "Studies 
of Sounding and Imaging Measurements from Geostationary Satellites." 
Significant progress has been made on most tasks during the reporting 
period. The prototype McIDAS equipment has continued to work well. 
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Introduction 

This third quarterly report covers work performed from 1 February 1973 

through 30 April 1973. In summary, progress by task is as follows: 

Task A . Investigations of Meteorological Data Processing 
Techniques 

Substantial progress has been made towards correcting ATS line 
start errors. This includes the development of an algorithm to 
detect the right hand earth limb in very noisy ATS images and the 
development of a processing technique which uses the limb detection 
algorithms to correct raw ATS data. Line start errors have been 
corrected on 23 ATS images, thus allowing relatively accurate wind 
determination from otherwise unusable data. 

Task B . Sun Glitter 

Mr. Jack Kornfield, the student who has almost completed a Ph. D. 
thesis on the "Determination of Sea Surface Wind and Stress from 
Sunglint Observed by a Geostationary Satellite," is no longer at the 
University. However, he is still working to complete the thesis. 

In order to meet the University thesis requirements, it must be 
completed by October. Therefore, it will be available as a final 
report at that time — whether completed as a thesis or not. 

Task D . Cloud Growth Rate 

In a master's thesis now completed, the effects of changing 
sun-satellite-cloud geometry on ATS-3 viewed cloud brightness have 
been documented. The brightness contrast analyses within clouds 
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were found to depend on scattering angles. For a large cloud mass 
different parts of the same clouds appeared to have different 
scattering properties; hence, different normalization corrections 
may be required. This suggests an important deficiency of present 
normalization procedures, which apply the same correction factor for 
all parts of the cloud. 

Also considered were observations of the scattering phase function 
and the effects of cloud breaks smaller than camera resolution. Work 
is underway on the development of a normalization procedure which will 
use a theoretical multiple scattering program to generate data on the 
brightness of clouds with various thicknesses for varying sun-cloud- 
satellite geometries. 

Task E . Comparative Studies in Satellite Stability 
Appendices 6 and 7 and Part VI to Das' technical report have been 
completed. Part VI of the report is attached at the end of this 
summary as an appendix. Complete coding of the model is proceeding 
on schedule. 

Task F . High Resolution Optics Study 

The seminar on mesoscale monitoring and prediction has been completed. 
Students have studied the causes and behavior of tornadoes, hurricanes 
heavy snow swaths, thunderstorms, hail, sleet, flash floods, clear air 
turbulence, fog, and lake and sea breezes. Key observables in the 
motion, thermal, and moisture fields were identified and means of 
communicating such information to the local Madison community explored 
Two days of briefings were held to transmit this information from the 
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students to SSEC physicists and engineers who will complete the study 
by investigating how the key observables can be measured; and, if 
time permits, by broadening the scope of the seminar subjects. We 
have reached two main conclusions so far: (1) SEOS will combine 

attributes necessary for mesoscale study which are not found in any 
other observing system, and (2) a high resolution vertical temperature 
and moisture sounding capability (in addition to imaging) is required 
to adequately treat mesoscale monitoring and prediction. 

Operation of the McIDAS system has been virtually trouble-free. Use of 
the system has been so intensive, the two primary problems over the last 
few weeks have been scheduling users and near saturation of the disk 
storage capacity. Certain program objectives and emphases have been shifted 
to reflect the demands of the new Task F, High Resolution Optics Study. 

Task Progress 

Task A . Investigations of Meteorological Data Processing Techniques 

During the past quarter, effort on this task has resulted in significant 
progress in the correction of ATS line start errors. Specific advances include 

1. Development of an algorithm to detect the right hand earth 
limb in very noisy ATS images; 

2. Development of a processing technique which uses the limb 
detection algorithms to correct raw ATS data; and 

3. Correction of line start errors in 23 ATS images, thus allowing 
relatively accurate wind determinations from otherwise unuseable 
data. 


Each of these items is explained in the following series of discussions. 
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Development of Right Hand Algorithm 

Development of a noise compensated threshold detection algorithm 
for locating an effective left hand earth edge has already been reported 
on in previous quarterly reports. The difference between the left hand 
and right hand edges of (some) ATS pictures is the amount of communication 
noise present. This noise on the right hand side is so large (for days 
203, 204, of 1969) that an entirely new algorithm was needed. 

Line plots of ATS data (day 204, 1969) near the right hand limb 
revealed large and frequent noise spikes both within the earth's disk 
and in space. Application of a simple filter (10 sample running mean) 
eliminated the spikes but did not remove lower frequency noise components. 
However, further filtering was not practical because the rise of the 
limb radiance profile would also be filtered out. Instead, the problem 
of falsely triggering on a low frequency noise burst was largely avoided 
by requiring the threshold to be exceeded at several points over a wide 
sample spacing simultaneously. In other words, the edge position was 
found at sample N if N was the largest sample number for which the ATS 
signal exceeded the specified threshold at N, N-l, N-2, N-5, N-10, N-20, 
and N-30. (A better, though more complicated procedure, would be to 
find the best fit position of a theoretical profile to the data.) As in 
the case of the left hand algorithm, the threshold is augmented by the 
mean space noise before testing against the data. 

Test runs with thresholds of 5, 10, 20, 30, and 40 were made to find 
the threshold which yielded the smallest line to line variance. The 30 
threshold was found to be the most stable. The lower thresholds appeared 
to be quite sensitive to noise; and the higher threshold appeared to respond 
to cloud cover variations. Unfortunately, it was not possible to check the 
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threshold detection results against landmark measurements (as was done 
with the left hand algorithm) since no suitable landmarks were illuminated 
at the same time the right hand limb was illuminated. The only further 
verification available at this time is in the consistency of wind fields 
derived from data corrected from the algorithm results. 

Development of a Data Processing Technique 

Given a measured position of the earth edge for each line of each 
ATS image in a sequence, there are a number of ways in which that 
information can be used to improve the accuracy of cloud motion vectors. 

The method we have chosen after some experimentation is one which maintains 
maximum compatibility with already existing data processing software. The 
basic operations involved can be summarized as follows: 

a) image navigation, 

b) calculation of navigation predictions of earth edge 
position as a function of line number for each image, 

c) threshold detection of earth edge position as a function 
of line number for each image, and 

d) write new digital tapes which are identical to raw tapes 
except for line shifts equal to the differences between 
results of b) and c). 

After d) is completed, images can be displayed and processed for wind 
sets using standard techniques. 

The first step (image navigation), although a standard procedure in ATS 
data processing, required special procedures and software to allow use of 
earth edge corrections of 'line start timing errors. 
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The present procedures are being employed in processing data on 
the WINDCO (prototype McIDAS) system. We plan to employ somewhat different 
and more efficient procedures on the advanced McIDAS system. These will 
be discussed in subsequent reports. 

Data Processing Activities 

In response to needs of Center scientists, corrected digital tapes 
have been written for ATS-3 data from day 204, 1969 and day 203, 1969. 

The corrected images are 800 line swaths encompassing the target area of 
the experiment. We now have 16 such corrected images for day 204, ten 
based on left hand limb detection and six based on right hand limb 
detection. In addition, seven corrected images have been written for day 
203, these based on left hand limb detection. Additional 203 data, as 
well as 202 data, will be processed during the next quarter. 

Wind sets have now been obtained for both right hand and left hand 
based corrections. They show east-west residuals which are comparable 
to, or less than, those obtained from raw data with low line start timing 
errors. This preliminary verification that the correction procedures do ■ 
result in significant data improvement will be made more quantitative 
by further tests planned for the next quarter. 
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Task B . Sun Glitter 

Sun glitter studies have followed two paths: (1) an investigation of the 

utility and technique of obtaining sun glitter measurements from a geostationary 
satellite images, and (2) a study of the problems of geometry (pointing angle, 
etc.) from a near earth orbiting satellite. 

Kornfield has almost completed a Ph.D. thesis on the first subject. He 
shows that one can indeed obtain excellent estimates of surface wind velocity 
(we already knew we could get good estimates of wind speed) from sun glitter 
measurements. There are two ways to .get the wind direction information. One 
can use the mean square slope in the cross wind and the downwind directions 
separately. The sunglitter intensity is a measure of the slope probability — 
but it is also affected by the intervening atmosphere, scatter by the white 
caps, scatter by the ocean, and also glitter arising from the sky dome. He 
shows how to account for all of these and is able to get wind direction for 
typical conditions, but not the very light winds or very strong winds. He has 
also found — and this may be more important — that the displacement of the most 
probable slope; i.e., the brightest part of the sun glitter, from the position 
the sun's image would occupy if the sea were calm, is directly proportional to 
the square of surface wind. Indeed, it appears to be directly proportional to 
the wind stress on the sea surface. This kind of information is difficult to 
extract from geostationary satellite images because they cover such a wide area. 
This data would be easy to extract from near earth orbiting satellites. (Kornfield 
started this work with support from another sponsor, but his support until recently 
has been under this contract. Although no longer at the University, Kornfield 
is continuing to write his thesis, which must by University rules be submittedly 
by October 1973). There is no specific progress to report on item (2) — the 



8 


geometry for the near earth satellite — except that the computer programming 
effort continues. 

Task D . Cloud Growth Rate 

A master's thesis outlining the requirements of a normalization procedure 
for ATS satellite picture has been completed. The draft Table of Contents 
and Introduction of the thesis is included below. In this work, using ATS-3 
data the effects of changing sun-satellite-cloud geometry on the cloud brightness 
have been clearly demonstrated. One of the primary deficiencies of previous 
brightness normalization procedures has been that they apply the same correction 
factor for all parts of the cloud. The brightness contrast analyses within 
clouds are found to be a function of scattering angles. For a large cloud mass 
different parts of the same clouds appear to have different scattering pro- 
perties and therefore may require different normalization corrections. 

This study also looked into the observations of the scattering phase 
function and the effect of possible breaks in 'the clouds which were smaller 
than the camera resolution. The results of this investigation in the form of 
the completed thesis will be included in the annual report. Work is presently 
proceeding on developing a normalization procedure which will use a theoretical 
multiple scattering program to generate data on the brightness of clouds with 
various thickness for varying sun- cloud-satellite geometries. 
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The purpose of this study is to look into the requirements of 
a normalization and standardization proceedure for satellite pictures. 
Meteorological satellites now provide pictures of the earth and its 
cloud cover on a routine basis. These satellites offer the coverage 
and resolution necessary to study and monitor mesoscale phenomena such 
as convection.. The satellites presently provide pretty pictures. 
However they offer the possibilities of providing considerablely 
more information on mesoscale phenomena than is presently obtainable. 

If one is to try to study mesocale phenomena from a satellite, one 
must first understand how light interacts with a cloud and what 
information about the clouds can be gotten from the light which the 


satellite receives. 


The reflected light from clouds depends to a large extent on 
four variables: 

a) Drop size distribution and phase state of water particles 

b) Number density of scattering particles in the cloud 

c) Cloud thickness 

d) Angular condition of the measurement system (the zenith angles 
of the sun and sensor and their relative azimuth angle) 

If one is to study some convective process, such as cloud growth, 
one wants only one variable to change at a time. If the cloud 
thickness changes, but the sun or the satellite also moves, one needs 
a method of removing the effect of the change in the measurement 
geometry. This is the normalization problem. Once the normalization 
problem has been satisfactorily solved, a considerable amount of 
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infopma^iQh tca-h^be gotten from satellite photographs . The volume 
flux through a cloud can accurately be determined along with the 
associated rainfall and energy release. This combined with the 
extensive coverage offered by the satellite provides possibilities 
for exploring the large scale energetics associated with convective 
systems o But first the normalization problem must be solved. 

There have been several approaches made at the normalization 
problem. The simplest approach has been to assume that clouds are 
perfect isotropic reflectors and obey Lambert's law. The intensity 
of the reflected light will vary as the cosine of the solar zenith 
angle for an unchanging cloud. There has been some evidence that 
this can be a valid assumption for very thick clouds. Martin and 
Suomi '(72) have shown that the tops of cumulonimbus clouds display 
a Lambertian behavior. Sikdar and Suomi- (72) also showed that thick 
clouds behave closely as Lambertian reflectors for small solar zenith 
angles (+30° to -30°) . 

However there is considerable evidence in the literature that 
other clouds do not behave as Lambertian reflectors. Bartman(67) 
measured the scattering off stratocumulus clouds and found anisotropic 
reflectance. Ruff, et. at, (68) measured the angular distribution of 


solar radiation reflected from clouds as determined from TIROS IV 


, radiometer measurements. They found clouds generally show an anisotropic 
reflectance pattern which varies with solar zenith angle. Erennan and 
Bandeen (70) measured the pattern of reflectance of solar radiation 
from stratocumulus clouds using an aircraft-torn medium resolution 
radiometer. Their results show clouds to be anisotropic with the 
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Sine© there is ample evidence that clouds are not isotropic 
reflectors, Sikula and Vonder Haar (72) . developed a normalization 
proceedure based on the empirical data of Brennan and Eandeen (70). 

This proceedure used empirical anisotropic factors which were tabulated 
for various solar zenith angles, satellite zenith angles, and the 
azmiuth angle between them. This proceedure is an improvement over 
the isotropic assumption, but it still has several drawbacks. One 


is that it needs continuous revision as more data becomes available, 


Another is that it makes no allowance for variations in the types and 
variations in the structure of clouds. For any given sun-cloud-satellite 
geometry, the same normalization factor is applied to all types of 
clouds, and to all parts of a single cloud, A cloud is a cloud, is a 
cloud; or is it? 

The purpose of this study will be to -look into the requirements 
of a letter normalization system and to provide some suggestions 
as to possible methods of improving the normalization prcceedures. 

In order to accomplish this, first the geometry of scattering and 
the definition of terms will be reviewed. Then the theoretical 


literature will be reviewed. The dependence of the single scattering 
phase function on the scatter's size distribution and shape will 
be briefly covered. For multiple scattering, the effects of the 
cloud's thickness, size distribution, and phase function shape 
will be covered. Then experimental observations of phase functions 
and of real clouds will be reviewed. In order to provide answers to 
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questions rb-isgd. by t.‘he“'theoretical and experimental literature 


review, I have made some observations. The question of does the 
contrast between the brightest parts and other parts of a cloud 
depend upon scattering geometry has not been answered by the 
experimental literature even though the theoretical literature, 
suggests that the contrast does depend upon scattering geometry. 
This question is important for normalization because one must know 
if the same normalization factor can be applied to all parts of a 
single cloud, or if different factors must be applied to different 
parts. I will try to answer this question and also the question 
of v/hat types of phase functions are found in clouds. Included in 
my analysis will be a discussion of possible effects, such as data 
biasing and the effect of possible inhomogeneties in the clouds 
which could cause misinterpretations of my results. 
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Task E . Comparative Studies in Satellite Stability 

Debugging and testing of the computer coded stability model continues. 

Part VI, "Boundary Layer Response and Asymptotic Stability" is attached as 
an appendix to this progress report. 

Task F . High Resolution Optics Study 

In February 1973, NASA requested SSEC to reprogram part of its activities 
toward a study of potential meteorological uses of the Synchronous Earth 
Observation Satellite (SEOS) . The basic purpose of the system would be to 
detect and predict hazards to life, property, or the quality of the environ- 
ment, and to promote proper exploitation and conservation of resources. At 
that time, a seminar on mesoscale monitoring and prediction has been completed under 
the direction of Professor Suomi. Topics of study were chosen with an eye 
toward possible SEOS applications, but the students were interested mainly in 
what observables should be observed and in applying them to improving local 
weather communication to the general public. 

In this seminar, the Causes and behavior of tornadoes, hurricanes, heavy 
snow swaths, thunderstorms, hail, sleet, flash floods, clear air turbulence, 
fog, and lake and sea breezes were summarized. On the basis of these summaries 
the key observables in the motion, thermal, and moisture fields were tabulated. 
Individual charts were made and combined to form a composite chart identifying 
common observables. The information about the key observables to be monitored ■ 
has been presented to SSEC project physicists and engineers in two days of 
briefings shortly before the end of the semester. The SSEC people will now 
consider how the measurements should be made and consider the impact of these 


requirements on possible sensing systems, especially SEOS. 
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As a result of post-briefing discussion, two primary conclusions are 
already evident: 

1. SEOS in uniquely important, and virtually indispensible , for monitor- 
ing and study of the meteorological mesoscale. This is because it is 
the only observing platform combining the attributes of large area 
coverage, high resolution, high signal to noise ratio, and continual 
temporal coverage. Other observing platforms fail in one or more of 
the above to satisfy observing requirements. 

2. SEOS will need a vertical temperature and moisture sounding capability. 
The combination of high resolution sounding and imaging is far more 
powerful than either alone since by knowing the temperature and moisture 
fields as well as the motion field one can get a handle on the dynamics 
of mesoscale activity as well as the kinematics. Knowing the forcing 
function permits inferences about future change, the key to prediction. 
Imaging alone tells us primarily what has happened, not what may happen. 

New Technology 

No items of new technology have been produced during this quarter. 

Program for next reporting interval 

Task A . Investigations of Meteorlogical Data Processing Techniques 

The right and left edge correction procedure will be applied to tapes for 
day 202, 1969. The finding that this convection procedure does result in 
significant data improvement will be made more quantitative through tests 
beyond the simple comparisons of east-west residuals with residuals from raw 
data of low line start timing errors. 
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Additional tests will be preformed to establish relationships between the 
level in the atmosphere where brightness values indicate detection (limb 
radiance profiles) and sun satellite geometry, clouds at the limb at the place 
of scan, and latitude. 

Task B . Sun Glitter 

This task has been curtailed to provide funds for the new Task F. Kornfield's 
completed Ph.D. thesis should be available by October 1973; in any case the 
work will be made available as a final report. 

Task D. Cloud Growth Rate 

Work will begin on a new brightness normalization procedure which will use 
a theoretical multiple scattering program to generate data on cloud brightness 
as a function of thickness and varying sun-cioud-satelli te geometries. 

Task E . Comparative Studies in Satellite Stability 

Debugging and testing of the computer coded model will continue. 

Task F . High Resolution Optics Study 

In addition to the summary charts prepared by the mesoscale monitoring 
seminar mentioned above (charts which stress the meteorological factors) , a 
systematic classification matrix is being assembled in order to optimize the 
specific technical requirements which must be met to successfully monitor the 
meteorological observables. Measurement Requirement sheets such as the sample 
attached will be filled out and the information transferred to punched cards 
and computer correlated and classified to develop a consistent set of technical 
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requirements which is complete for the phenomena we have investigated. 

If time permits, the study will be broadened to other mesoscale phenomena 
and potential uses of SEOS in order to make a more complete survey. For the 
present, we are concentrating on analyzing and documenting the unique and 
significant contributions of the seminar. A full description of the study 
results will be completed before the next quarterly report. 

Conclusions and Recommendations 


None at this time. 
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Part VI 

" Boundary Layer Response and Asymptotic Stability " 

1 . Introduction : 

In the following, the term "Boundary Layer Response" refers to the motion 
of a flexible satellite near the time t = 0 , and is not to be confused with the 
response of the atmosphere. 

In the earlier phases of this work, it was seen that the response of the 
system near the time t = 0 was found to be given essentially by the outer 
boundary layer equations (3.70) and (3.71). These are rewritten here as follows 

A(w b ,0)£(oj b ,t) + /e B(w b , 0 )£(w b >t) + £(w. b > t ) = 0 (3.70; 6.1) 

[AVyJ (w B ,0)io (t) + 2[ A7y ]((U B ,0)h)^(x) 

+ _7[Ay+Y.] (a) B ,0 )oj^ = d^(m B ,0) (3.71; 6.2) 

It was also seen that the long time response of the system is basically 
governed by Eqs. (5.58) and (5.113). These equations are 

. Bx + Bx+Bx = 0 (5.58; 6.3) 

1 2 3~ 

and 

. A q + A q + A q = A (5.113; 6.4) 

10 20 o 30 o —40 

Equations (6.1) and (6.2) will now be considered in detail for the 
boundary layer response. Then Eqs. (6.1) through (6.4) will be examined for 
stability of the satellite motion. The notation in these equations were 
defined before. 

2. The Boundary Layer Equations 

The functional matrices A and B of Eqs. (6.1) and (6.2) are defined in 
Eqs. (3.6) and (3.7). In Eqs. (6.1) and (6.2) these matrices are to be 
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evaluated at t = 0. This makes these two linear second order equations 
with constant coefficients. Hence _jr(x) is given by 


£(l.) = ^(0) + ^(O) 


(6.5) 


where <j> and <f> are fundamental matrices of Eq. (6.1). Thus y(x) 
3 4 

undergoes damped oscillations. For the two bodies A and B, Eq. (6.5) 
becomes 


—A3 ‘ *A3 “a 3 (0) + V i?A3<°> 


( 6 . 6 ) 


»B3 - *B3i!B3 (0) + V, V 0) 


(6.7) 


Substituting the forms of and ^ from Eqs. (4.21) and (4.23) 

into Eqs. (6.6) and (6.7), respectively, we obtain 


S A1^B + S A1® + 2S A1 W B + S A2 + S A3 + S AsVa2 


- S A5t 6 A2% ] + Wa ~ 6 A [I A W B ] 1 + i S Al (W B 


+ V B ) + S A2“B + S A4 + S A5 W B 6 A2 l0 B " WW^B 


+ VbVb " 6 A [I A^B ]W B + S A2% 


+ t S A5 9 V + 6 A 9I A® ] + ( S A5 06 A 2 V 

+S A5 iVl«A2 ^ + VaV + 


+ < s a 5 [ “i^ I6 A2V + 6 A [ V !I aV 


+ < S, a> + d.ULl.ur, \ - S.u. 

| A5 B A2— B A B A— B ( A— A 


+A V 0) + VA 3 <0 > 


( 6 . 8 ) 
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and 


S B1-B + ' S B2^ + S BA + S b4 + S B5“B S B2^B 


+6 bVb% " Vb = ^B3~B3 (0) + ^B4^-B3 (0) 


(6.9) 


From Eq. (6.9), _6 is solved in terms of oj , and we obtain 


6 = expt-s^s B4 t] 




eXp[S B2 S B4 T] tE - (T) 


*B3^B3 (0) + *bA 3 ( ° )] 


( 6 . 10 ) 


Eliminating Q_, J3 and ^ from Eq. (6.8) by using Eqs. (6.9) and (6.10), the 
boundary layer equation for u> , corresponding to Eq. (4.6a) is obtained as 

— D 


fi l^B + B 2 ~B + B A + \) I B 5 (T ' t) ^B— B 


! 


II 


■ S lA - S B3SB )<t)d h + B 6~3 + ®7— B - 6 A±A 


+ E ’ite B , i o B ,JJ | U A ,u B ) = B s (t) 


( 6 . 11 ) 


where. 


B (t) = -B 4 F x (t) 


0 - 


+ (S A2 +S A3 )S B2^B3 (t) 


4 3 (0) 


+ < S A2 +S A3 )S B2V (t> ]2B3'°> 

— B A3j [*B3 (T) 433 C °> + V°*B3< 0) 


+ *A ^AS™ + +A3 (t) i!A3 (0> : 


( 6 . 12 ) 


As the body B is nominally static, the terms involving the products of 

03 or w with 1L (t) have been neglected in deriving Eq. (6.11). 

— B B B 3 
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3. Boundary Layer Solution for to 


£1 


■B 


The maximal probabilistic boundary layer solution for to is obtained in 

— D 

a way similar to that used for the asymptotic case. Corresponding to Eqs . 
(5.44) and (5.45), the governing equations in this case become 


and 


T„ T T 

Vl - B 2— 1 + B jh " £ ' 


B l-B + B 2-B + B A ' \ 


g A +g A 
& 1-1 6 2 — 1 


(6.13) 




(T)dT ' 


+»(t) - B s (t) - e' 


g,(t t)A (T)dT + g A + g A - f 

3 1 1 4 1 5 1 — 


(6.14) 


Comparing Eqs. (6.13) and (6.14) with Eqs. (5.44) and (5.45) we see that 

k 

the term B F 0 of Eq. (5.45) has been replaced by -B (t) in Eq. (6.14) 

4 1 8 

without any other change. Hence all the expressions for to , derived for the 

... B 

asymptotic case becomes applicable to the boundary layer case if -B (t) is 

8 

* 

substituted for B. F ,0 . 

4 1— 



The governing equation for the boundary layer solution for is given 

by Eq. (6.2). As the coefficients of and to^ are evaluated at 

t = 0, the equation is one with constant coefficients. The initial conditions 
are lumped on Eq. (6.1) and so has zero initial conditions. Hence the 

solution for the bodies A and B are given respectively by the Duhamel’s 

integrals as 


t 

r 


— A1 



(t-T)d(o) B> 0)dx 


(6.15) 
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and t 

r 




-Bl 


<f>g 5 (t-x) d.(ojg ,0) dr . 


(6.16) 


In the above two equations, <f> and <f> are the appropriate fundamental 

AO ij o 

matrices of the homogeneous equations of Eq. (6.2) for the bodies A and B, 
respectively. 


5 . Stability Criteria 

The basic requirement for a stable satellite motion is that the angular 

velocities w and u> should approach zero values as t -*■ 00 . The 
B , 1 B , ^ A 

components x. , (i=0, 1,2, . . . ) , of w are given by Eqs. (5.47) and (5.51) 

through (5.56). It is seen that the forms of the solutions for x. are 
governed by the roots of the homogeneous equation (5.58). For convenience, 

Eq. (5.58) is repeated here as 

B,x/ + B x. + B_x. = 0 (5.58; 6.17) 

l—i 2—i 


For asymptotic stability, the roots of this equation must have negative real 
parts. The characteristic equation of Eq. (6.17) is 


p b + tr-[B 1 B 2 ]p D + jtr-lBjBg] + DetlB^tr- [B^] p H j 
+ Det[B"|]|Det[B 2 J + = rs .rlh ,sJ B 2 , tk| p ' + j^t- [B^B,] tr- [3,3^ 

‘ J ' 


+Det* [B 


3jB 2 ]tr- [B 2 B 3 ]ip 2 + |Det.[B 1 1 B 3 ]tr-[B 1 3 B 2 ](p 


-1 


+Det[B 1 B 3 ] = 0 


(6.18) 


In Eq. (6.18), e and e are the permutation symbols. 

1 j k. ITS t 

The Routh series for Eq. (6.18) is now calculated and hence we obtain 
the following inequalities [1] to be satisfied for stability of the linear 
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core of Eq. (4.46) . 

-1 


73 


tr. [B^] > 0 

tr- [B J b 2 ] -jtr. [B 3 J + Det • [B ^ ] tr • [B^B 1 ] ' 


(6.19) 


-Det [B j ] | Det[B 2 ] + 


e e . B, . B , ,B 0 . / 

rst ljk d ,n i ,sj 2,tk| 


\> o 


( 6 . 20 ) 


DetfB^ *|tr- [B^B 2 ]| |tr*[BjB 3 ] + Det [B |b 2 ] tr • [BgB^j |oet[B 2 ] 


+e rst 6 ij k B 3 , ri B 1 , s j ® 3 , tk 


6 rs t S i j k^ 3 , ri^ 1 , s j ^ 3 , tk 


+Det-[B|B 2 ]tr.[B 2 B 3 ]^ + ^tr-fB^j ^ 




Det (B j ] |oet [B 2 ] + 

■ |tr-[B^B 2 ]| | Det • [B J b 3 ] tr • [B^B ^ ] 


mj! ! Det [B ^ 


~1 


B 3 ]tr- [B 3 B 2 ]>> 0 (6.21) 


Rewriting Eq. (6.18) as 

,6 j_ „5 -j. „ „4 j_ 0 r'. 3 


P + + a 2 p 4 + a^> i + a^p- + a 5 P + a 6 = 0 


(6.22) 


the fourth in equality is obtained as 


k=(aaaa +aaa - a 2 a - a 2 a 3 - a 2 + 2a a a 
1234 235 34 14 5 145 

- a i a 2 a 5 " a l a 2 a 6 “ a i a 3 a 6 ) > ° 


(6.23) 


The fifth inequality is given by 

[k (a i a 2 a 5 - a 3 a 5 - a^) + 2a 1 a 3 a 6 b i b 3 - a^b^ - > 0 


where 


b l = a i a 2 " a 3 


(6.24) 


(6.25) 


and 


b = a a - a , 
3 14 5 


(6.26) 


c 
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The last inequality to be satisfied is 


?A 


a g = Det- I b 1 B 3 ] > 0 


(6.27) 


The above inequalities will all be obviously satisfied in a well-designed 
satellite. These also yield a comparison technique for two satellites. The 
comparison is to be made on marginal stability. The design for which a^ and 
a are larger has a larger marginal stability and hence is the more desirable 
one. 

We now consider Eq . (4.25) to examine the conditions for stability of the 
relative rotation, 8. Equation (4.25) is repeated here for ease of reference. 



_i 

* 

/•t 

~ - 1 

= exp 

_S B2 S B4 t 

_9 + 

exp 

S B2 S B4 ( ‘ T " t ^ 

i 


c 

3 



g(t)dx 


(4.25; 6.28) 


where g(x) is defined by Eq. (4.26). Evidently, the primary requirement for 
stability of _6 is that the eigenvalues of the matrix [-S must have 

non-positive real parts. There is another interesting requirement on these 
eigenvalues. This arises due to the appearance of g(x) in the integrand in 
Eq, (6.28). This requirement is obtained in the following discussion. 

Let "Vjj - h 2 and -V> 3 be the real parts of the eigenvalues of [-S B * S g 4 ]> 
such that 


I V ! I > l^ 2 l > Ih 3 l 


(6.29) 


Let -a ,-a ,-o be the real parts of the eigenvalues of Eq . (6.17), such that 
12 3 


Iffjl > l a 2 l > l° 3 l 


(6.30) 


Then from Eq . (6.28), we have 

r t 

~ y 3 t I V _y 3 t _a 3 X 

0 [ < 8e + I a * e • e • e dx 


o 



25 


8 


"V “l - ( °3 +M 3- p l )t: 
0 < 6e + e 


^r°3 


(6.31) 


Hence for J0 to be asymptotically stable,, we must have 


°i 4 1=1,2, 3. 


° 3 — - (y l" p 3 ) * 


(6.32a) 


(6.32b) 


The variables 0 and introduced in Eq. (6.31) are suitable finite con- 

stants . 

We now consider Eq. (5.60) for the stability of the component x (t) of 

o 

<jdu(t). For convenience, Eq. (5.60) is rewritten below. 

£> 


f ( 

X Q (t) = ^(t)^ + <f> 2 (t)b 2Q - <f» 2 (t-x)B J F^t: 

J ' 


-j B 5 (S-t) [S b1 x q (0) + S B3 x q ( 0) ]ds| dx - J <f> 2 (t-x)w Q (x) dx. (5.60; 6.33) 

Since the first two terms and the last integral on the right hand side of 
Eq. (6.33) tend to zero as t -*■ the truncated form of Eq. (6.33) given 
below is considered. 


,>) = 'J <!> 2 (t-x)B 4 j Fl (x: 


- V s 


- T) [ S Bl^o (0)+S B^o (0) dS | dT - 


Going through the derivation of Eq. (6.33) it can be seen that the above 
equation is generated from 


^o (t) 


L (t-x)B 0 (t) dx 
z if- 


(6.34) 



Hence using Eq. (6.31), Eq. (6.34) is converted to 

„t 


, , . -° 3 t <V 

J2i 0 (t) | 1 a 2 -| e 


o 


"V (y r H 3 -a 3 )T 
e +a 3 e 


or 


a 2 a 3 


< A_. e ( VV°P< 

~° '-a r u 3 (a 1 +y 1 -a 3 -y 3 ) 


dr 

(o 1 +V 1 -P 3 -20 3 )t 

1 

(6.35) 


From Eq. (6.35), the necessary stability conditions are 


°1 -°3 +y 3 


°1 1 2 a 3 + y 3 -^! 


(6.36) 

(6.37) 


Inequalities (6.32b), (6.36) and (6.37) can be combined to obtain the following 
system of inequalities. 


° i .i..( a 3 + M 3 ) 5 <? 3 > V, (6.38) 

or 

a ! 1 ( 2 cr 3 +y 3 _ y i) > ]i 1 > 0 3 > (b 1 -y 3 ) (6.39) 

t 

If the above inequalities are satisfied, both j) and Xq tend to zero 

r 

asymptotically. It is to be noted that 2L 0 ( t ^> given by Eq . (6.33), is also 
the solution of the linear core of the deterministic equation (4.46a). Using 
this fact, we now prove that the non-linear equation (4.46a) is stable at the 
origin. 

/ . 

The non-linear term in Eq. (4.46a) is e f (or, ,u> ,u. ,u ) . Comparing Eqs. 
(4.46a) and (4.35), it can be seen that e _f consists of quadratic and higher 


order terms of its arguments. Since 



we have 


Lim 

L. 

0) -> 

0 

-B 


ID -> 

0 

-B 


V = 1 

fi 


ef (v) = 0 


v 




(6.40) 


-LJLOU 


Then from the Theorem 1, of Struble [2], and the non-linearity condil 
(6.40), it follows that Eq. (4.46a) is stable at v(t) = 0. 

We now analyze the behaviour of the stochastic differential equations 
(5.44) and (5.45). These non-linear equations are 


T" T* T ' 

B,A, - B^A + BY = e 
1—1 2 3-1 


8 

t 


(5.44; 6.41) 


B l“ B + ®2— B + B 3^B “ \ {\ B 5 (T - t)(S BlV S B^B )(T)dT | 

t 


+w(t) + B^F^e = e 


11 


8 3 (t,T)A 1 (0)dT + 85 ^ “ JL 


(5.45; 6.42) 


0 

It can be proved that the eigenvalues of Eq. (6.17) and those of the left hand 
side of Eq. (6.41) are equal in magnitude but opposite in sign. Thus A_^ is 
unbounded. This result is expected as _A , which is a measure of the indeter- 
minacy of the value of should grow larger with time. An interesting 

point to be noted is that the larger the marginal stability of co^, the larger 
is the marginal instability of A_ . Equation (6.42) is now considered to 
resolve this apparent paradox. 
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As it stands, Eq. (6.42) gives unbounded solutions for w . This can 

— B 

be seen from the following. Let only the term t> e considered on the 

right hand side of Eq. (6.42). The solution corresponding to this term is 
given by 


x(t) = e 


Since 

g (t) = 
5 


<f> (t-T)g 5 (t)A 1 (x)dT 




/ 3 1 d 9 l\ T " 

B 6di + B 7 


B Bl3u dt 3u / 



L \ — B -B / J 


* 3f 
6 ,U 


A A 3u. 
—A 




g c (t) behaves like e 

O 

Hence d 

, , %l , ,( -V+V < 0 f°3> T 
| x(t) I 1 ke 


dx 


(6.43) 


or 


2£(t) _< k^ 


2(Oj-a )t 


(6.44) 


where k and k-^ are suitable finite constants. Thus from Eqs. (6.30) and 
(6.44), 

|x(t) | -+• OT as t ->■ 


To have bounded solutions for , we must have either 


°1 = ‘ S " ° 


(6.45) 


or 


°A - U B ■ ° 


(6.46) 


or 


Det 


" 31 

J 

3f ' 


" 3 1~ 



d 


= Det. 


_ 3 — B 

dt 

B _ 

3 u . 
— A_ 


(6.47) 



12 


The condition (6.45) ensures purely oscillatory motion of the satellite and 
is too restrictive. The condition (6.46) implies deterministic forcing and 
control torques. This also is quite difficult to achieve. The implications 
of the third condition (6.47) will be analyzed now. 

Expanding the non-linear terms of Eq. (5.45), it is seen that _f does 


not contain u„ and u.. Hence 
— B —A 


’ 3f ' 

:• ■ji — - ° 


(6.48) 


— - 0 
3 ^b 


(6.49) 


We also have 


-r — = - 2S. .0) + S.,u) 6 Ao - S . (6 w ) + 6m I A 

3u A1 B A5 B A2 A5 A2— B ABA 

— B 

-°A (I A^B ) S B2'i + S A5 (S B2 S B1-B ) + <S B2 S B3% > 


+<S B2 S B5"b 6 B2“ b ) + ( S B2 S B“ B I B2B ) " 


+ < S B2 S B^ S A2 S B2 S B + S A5 f < 6 A2 S B2 E B + < S A 2 S m S bA> 


-1 


-1 * 


-1 * 


+ < 5 A2 S B2 S B5"b 5 B 2 V + < 6 A 2 S B2 V B W ' A.Vb^ 


+ (6 A2 S B2 S B^ S B2^B + (S B2 S bA> + ‘V B A> 


+ ^ S B2 S B5° J B <S B2-B' ) + ^ ^ S B2 6 B-B^ 


^ ^JL W j Q 

+ < S B2 S B 4> Wb + S A (I A S B2 S bA ) + (I A S B 2 S B ft> 



13 




+ ( ‘ I A S B2 S B5 a) B 6 B2-B' > + ( ‘ I A S B2 5 B^B I B^B' > ( ‘ I A S B2 6 B-B' ) 




+ ( I A S B2 S b4> 


-1 * 


S B2 6 B - S A5< 6 A2“B® )S B2 6 B " S A5 < V> 6 A2 S B2 6 B 


* 


-WUaV’Ob - 6 1 ( "&® I A S B2 5 B 


(6.50) 


From the above equation, it is seen that Eq. (6.47) can be satisfied if 


Det 


-1 * 
S B2 6 B 


= 0 


(6.51) 


Since Eq. (6.51) can never be satisfied, it seem that the random torques tend 

to increase the variances of to indefinitely though the mean values of to 

— B ‘ D 

decrease asymptotically. 

6 . Single Composite Body Analysis 
a) Equation of motion 

We now modify the previous analysis for a satellite which consists of the 
body B only. This satellite may either be a spin stabilized or a three-axes 
stabilized one. All equations through Eq. (4.3) remain valid. Since there 
are no contact forces and torques for a single body, equations corresponding 
to Eqs. (4.4) and (4.5) become 


A = [6 B1 ]a) B + VS B2 ]m B 


P = u B - [I B ]m B + [h B ]co B - a> B [I B ]jS B 


(6.52) 

(6.53) 


Then the equation for to , corresponding to Eq . (4.23) is 

is 


S Bl— B + S B 3— B + S B5 W B [6 B2 ] ^B 


^WV^B ~ [6 B ] ^B = ° 


(6.54) 
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b) Deterministic solutions for id. and u„ 

It is known that in the torque-free rigid problem, a constant value of 
w is a solution. This constant solution vector is the nominal angular 

— D 

velocity vector of the satellite. With this information, a perturbation 
series is now assumed for to . Let 

-D 


h = e x + (e ; ) + ( e ' ) 3 + 


(6.55) 


and 


“ B = ^0 + * 


where t ^ ie const ant nominal angular velocity vector of the satellite, 

and 


£ = lu^Cx ) at t = 0 

D n 


(6.57) 


Let u„(x„+h) be expanded in a Taylor’s series, such that 

— D 0 — 


A 

6 b J1 b = e w x +. (e' ) 2 w 2 + (e' ) 3 w 3 + . . . 


(6.58) 


Substituting Eqs. (6.55), (6.56) and (6.58) into Eq. (6.54) and separating the 

I Q i 

coefficients of (e ) and e , we obtain the following equation: 


x + 
Bl— l 


S B3 +S B5 X 0 6 B5 S B5^ 6 B2-0' )+6 B X 0 I B 6 B < ' I B^0' ) 




=W, T 

—1 e 


S B^ 0 +S B5 X 0 5 B2 X 0 +6 B X o I B X 0 


(6.59) 


Equation (6.59) is rewritten as 


S Bli + (S B3 +S B6^1 * ^1 - F S 


_ 1 

e“B3^0 e 


S B7 ( -0 ) 


(6.60) 


/ O 

The coefficients of (c ) generate the equation 



32 
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s bA + ° ^ - 

From the coefficients of (e*) 3 , 

S B7 ( Sl ) ] Si 
we obtain 


S Bl-3 + ^ S B3 +S B6^— 3 “ — 3 

S B7&1> 

1 

1^' " 

1 

S B7^ X 2' ) 


(6.61) 


(6.62) 


Equations from higher powers of e can be obtained similarly. Equations 
thus obtained are linear with constant coefficients and will be solved as 
shown in the earlier parts of this work. 

Since the control torque functions w^ appear linearly in Eqs. (6.60), 
(6.61) etc., w^ can be shown to be a "band-bang" control by the method 
shown earlier. Thus the controls for a single body satellite are seen to be 
much simpler than those of the dual-spin configuration. 

c) Random solutions for w and u 

— X> — D 

We now follow the method outlined in Part 5 of this work to obtain the 
statistical parameters of the random torque responses. The functional, J, 
to be minimized in this case is given by 


J = 


Vo) - V°) 


t -ir 

R 

0 


V°)~- B (°) 


+ 


"(r - - 

T -1 

_ 

j z 1 (t)-z 1 (t) 

Rj (t) 

z 1 (t)-z' 1 (t) 

I *- -iT ; 

i r 

1 


B (t) Ji B (t) 


V* 


U B (t) 


u B (t)-u B (t) 


+2A, |S„. V S B A +S B5 U B 6 B2^B + VbVb 


- 6 b% 


+ A „ 


-^B~— E 


dt + 2A 4 (z^ 1 -h 1 ) 


(6.63) 


For J to be minimum, the following set of equations are to be solved simulta- 


neously : 
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S B1-B + S B3^B + S B5VB2^B + 


z^t) = li^WgCO) ,t 


6 #b = 0 


UgCt) = u B (t) + U B 


r *f 

M 


(t) 


T 

,(°) -V°> + r o[ s biAj«» + \] 


u B (t) = u B (t) + U B (t)X 2 


(6.54) 

(6.64) 

(6.65) 

( 6 . 66 ) 

(6.67) 


Zj(t) = z_j(t) + Rj (t) A_ 4 

a 


s bAi 


s„„ + 


_0)„6_ _0J +6„u> 


B 3 8u> b ^ B5 B B2— B B B B-Bj 
9 H B \T 


X 

-1 


X x (T) 


= 0 


( 6 . 68 ) 


(6.69) 

(6.69a) 


These equations are reduced to 


S B f aJ -B + S B^B + S B5 W B2^B + Vb^B 


= 6, 


H B <o+u B («‘) x l( t) 


and 


T • 

S B1— l - 


S B 3 +S B 5 W B 6 B 2- S B 5 ^ 6 B2-B ^ +6 B^B 1 B 6 B ( 1 B-B ^ 


-1 


(6.70) 


= 0 (6.71) 


To solve Eqs . (6.70) and (6.71) the following series are used: 


= x Q + e' -x' + (e 'f ^2 + ( £# ) V 3 + • • • (6.72) 

Xj = Aio + £ 'A n + ( £ ') 2 A 12 + (e ' )3 Ai 3 + • • • (6.73) 

* 

6^ = e'wj + (e'£w + (e') 3 w 3 + . . . 


(6.74) 



Let U be defined by 


e'u=6 B U B 6 B (6.75) 

Then from Eqs . (6.70) and (6.71) the following set of equations are obtained: 


T • T 


S B 1—10 <S B3 +S B 6> -10 0 





(6.76) 

S lA + <S B3 +S B6^'l ' - 

^S x -iT 
e B3^0 e 

5 B7 (^ 0 )] x 0 + 

u -h„ 

(6.77) 

S BlAll - < S B3 +S B6>\> ■ [ ! 


X 

~10 



(6.78) 

S B1~2 + ^ S B3 +S B6^— 2 “ — 2 

[ S B7 ( i 

■1>] £l + 

U ill 


(6.79) 

S B1-12 " ^ S B3 +S B6^--12 = 


ilo + 

s B 7 (—i^J in 


(6.80) 

S B1^3 + (S B3 +S B6^3 = ^3 ~ 

[ s b ■/"- 

—2 " 

2^| — 1 

+ w 12 

(6.81) 


Similarly, more equations of this sequence will be obtained. It is evident 

that the random solutions coincide with the deterministic solutions if 

IL, = 0. These equations also will be solved by the method outlined in Part 
a 

5. 

From this, the mean and the variances of the pointing error of any rigid 
body will be obtained as shown in Part 5. 
d) Stability criteria 

As in the case of the dual-spin configuration, the random parameter X_ 
is unbounded, if the mean values of the response are maintained stable. This 
also implies that the variances increase indefinitely. We now obtain the 
criteria for stability of the mean values of cm, . 

— D 



18 


From Eqs . (6.77), (6.79) and (6.81), it is clear that the asymptotic 
stability is governed by the roots of 

S Bll + (S B3 +S B6^’ ' 0 <6 ' 82) 

For a stable solution, the roots of Eq. (6.82) must have negative real parts. 
The characteristic equation of Eq . (6.82) is 

-M • P 3 + Det -[ S m]: tr [ S B! (S B3 +S B6 ) ] P 2 

+Cet ' [ S B3 +S B 6 ] [‘W^B'Jp + Det ' [ S B3 +S B 6 ] “ ° < 6 - 83) 

Using the Routh series of Eq. (6.83), the following inequalities are obtained 
as the required criteria: 


Det - U > 0 

(6.84) 

Det '[ S B3 +S B 6 ] * ° 

(6.85) 

“•[ S b! (S B3 +S B6>] > 0 

(6.86) 

“•[ S B 1 1 < S B3 +S B6 ) ]- tr '[ <S B3 +S B6 ) '‘ S Bl] > 1 

(6.87) 


7. Performance Criteria of Satellites 

Based on the preceding analysis we now propose several criteria for 
evaluating the performance of satellites of different designs. It has been 
shown that for each satellite there is an associated outer boundary layer. 
The duration of this boundary layer depends on the structural properties and 
the nominal angular velocities of the satellites. This duration is /e”, 
where Eq. (3.5) defines e as 
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e = the largest eigenvalue of 




(3.5) 


It was shown that the boundary layer is the region of turbulent motion of the 
satellite caused by its flexible elements. Hence our first desirable criterion 
is a low value of e. 

We now consider the motions of satellites beyond the boundary layer. 

Clearly, all useful information can be obtained only in this quiescent zone. 

It is assumed that all designs to be compared are stable, so that the quiscent 

oscillations will eventually die out in each case. Since the rate of decay 

may be different for different design's, we propose that a larger marginal 

stability of motion be taken as the second criterion for evaluation. 

/ 

Our third critical parameter is e , which is the maximum value of external 

/ 

torques acting on a satellite. This parameter, e governs the amplitudes of 
vibrations which should ideally be absent. 

The fourth and most important criterion that we are proposing is the average 
pointing error for a sensor in the interval between two consecutive measurements 
of the angular positions of the satellite. Let [0,T] be this interval. Then 


the average pointing errors or are defined by 


E Ai " T 


4i dt - 2o 


( 6 . 88 ) 


and 


E Bi ^ t 


4i dt 


(6.89) 


where ip,. and ip . are given by Eqs. (5.130) and (5.127), respectively. 


r Ai r Bi 

In the above the norm 


v 


of a vector v^, i = 1 - n, is given by 


n 

|v|l = ^ V 

i=l 


(6.90) 



8. Conclusions 


We now have come to the end of our analysis. In this analysis, the method 
of obtaining the motion of the sensors aboard a satellite is obtained by per- 
turbation techniques. The most probable motion is obtained after considering 
random torques, initial conditions and measurement errors. Several stability 
criteria have also been established. Finally, several performance evaluation 
parameters have been proposed. It is hoped that these parameters will be a 
useful tool in comparing the pointing accuracy of different satellites. 
Numerical results will be presented in the following parts of this work. 



